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Abstract 

We calculate analytically the current statistics for a short diffusive wire be- 
tween the normal reservoir and a short superconductor-normal metal-superconduc- 
tor (SNS) junction, at arbitrary applied voltages and temperatures. The cumulant- 
generating function oscillates with the phase difference (|) across the junction, 
approaching the normal-state value at (|) = 7t. At T = and at the applied voltage 
much smaller than the proximity gap A,],, the current noise P; doubles and the 
third current cumulant C3 is 4 times larger with respect to their values in the 
normal state; at eV ^ A,|, they acquire large excess components. At the gap edge, 
eV = A^, the effective transferred charge defined through dPj/dl and dPj/dV 
approaches Oe and 3e, respectively, which makes doubtful the interpretation of 
these quantities as physical elementary transferred charge. At T 7^ 0, C3 shows 
a non-monotonous voltage dependence with a dip near eV = A,^. 

During last years, statistics of quantum and thermal fluctuations of the elec- 
tric current in mesoscopic systems attracted rapidly growing attention. It was 
established both experimentally and theoretically that the fluctuation properties 
of mesoscopic conductors provide important information about correlations 
and statistics of charge carriers which are not accessible through conductance 
measurements. A powerful theoretical approach to the fluctuation problem has 
been developed by using the concept of full counting statistics (FCS), i.e., statis- 
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tics of number of particles transferred through the conductor. The concept of 
FCS, which first appeared in quantum optics, was extended to normal electron 
systems in Ref. [1], and then applied to superconducting structures in Ref. [2]. 

The basic problem of the FCS is to calculate a probability PtgiN) for A'^ 
particles to pass a system during an observation time to- Equivalently, one can 
find a cumulant generating function (CGF) S{x), 

exp[-5(x)] = £^P,o(iV)exp(Wx), (1) 

which determines the irreducible current correlation functions (cumulants) C„, 

Cn = -{d/idxrS{x\^,. (2) 

The first two cumulants, Ci =N = Lw^^^fo (^) and C2 = {N- NY, correspond 
to the average current / = {e/to)C\ and noise power Pj = {2e^/to)C2. Intensive 
studies of the current noise have led to a number of interesting results concerning 
statistical correlations in the current transport (for a review, see Ref. [3]), and an 
effective charge qejf transferred during an elementary transport event. The third 

cumulant = {N — has recently attracted a special interest as the lowest- 
order correlation function which is not masked by equilibrium fluctuations [4]. 
First measurements of Ci{V) in the tunnel junction [5] have revealed its high 
sensitivity to an electromagnetic environment [6]. 

In normal metal/superconducting (NS) hybrid structures, the basic mecha- 
nism of charge transport at subgap energies, £" < A, is due to Andreev reflection 
of quasiparticles at the NS boundary [7], i.e., conversion of electrons incident 
from the normal metal to retroreflected holes, accompanied by escape of Cooper 
pairs into the superconductor. During an elementary Andreev reflection event, 
the effective charge transferred through the NS interface is twice the electron 
charge, qeff = 2e. This charge doubling strongly affects the current statistics in 
NS junctions; in particular, it is known as the reason for the doubling, compared 
to normal junctions, of a zero-bias shot noise [2, 8]. At finite bias, the effec- 
tive charge becomes dependent on applied voltage [9, 10], due to variations 
of the size of the proximity region near the NS boundary, where the quantum 
coherence holds between the electrons and retroreflected holes. 

In the Andreev interferometers (see Fig. 1), the phase relations between 
the electron and hole wavefunctions in the normal wire can be controlled by 
magnetic flux enclosed by a superconducting loop, which results in a period- 
ical dependence of the transport characteristics of the interferometer on the 
superconducting phase difference (|) across the SNS junction. First, the oscilla- 
tions of the conductance have been investigated experimentally (see a review 
in Ref. [11]) and theoretically [12], and, recently, the oscillations of the current 
noise were reported [10]. 

Motivated by growing interest in higher correlation functions, we discuss 
in this paper the full statistics of charge transport in Andreev interferometers. 
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Figure 1. A model of Andreev interferometer: diffusive wire of the length L connects a normal 
reservoir (N) and short SNS junction of the length d\ magnetic flux <5 threads a superconducting 
loop (S) of the interferometer. 

We will adopt several simplifying assumptions, which enables us to present an 
analytical solution for the CGF and thus to clearly demonstrate essential features 
of coherent effects in the current statistics in NS structures. Our approach is 
based on the extended Keldysh-Green technique [13] (see also [14]) in which 
the CGF is determined by the equation 

{-ie/tQ)dS/dx = /(X), /(X) = ^ / dETrxJ, Xk = o,x,. (3) 

The Pauli matrices a (x) operate in the Nambu (Keldysh) space. The counting 
current /(x) is to be found from the quantum kinetic equations [15] for the 
4x4 matrix Keldysh-Green function G in the mesoscopic normal region of the 
interferometer confined between the reservoirs, 

(5n [c>-;£, G] = ih-D dl, I = ONGdG, G^ = l. (4) 

In this equation, -D is the diffusion coefficient, 3 denotes spatial derivative, and 
On is the normal conductivity per unit length. The counting field x is introduced 
via a modified boundary condition involving the gauge transformation of the 
local-equilibrium function Gr, e.g., in the right (/?) normal reservoir, 

Gr(x) =exp(/xx^/2)GRexp(-/xT^/2). (5) 

A brief overview of this technique in the particular case of normal structures 
can be found in the Appendix. 

For a multi-terminal structure in Fig. 1, the solution of Eq. (4) has to be 
found separately in each arm of the interferometer, with the matching condition 
following from the Kirchhoff 's rule for partial counting currents at the node 
[16]. The problem simplifies when the junction length d is much smaller than 
the length L of the interferometer wire (or, more precisely, the wire resistance 
dominates the net interferometer resistance). In this case, the interferometer 
wire weakly affects the spectrum of the junction [17], which thus can be consid- 
ered as an effective left (L) reservoir. Correspondingly, the function Gl which 
imposes the boundary condition to Eq. (4) at the junction node, is to be con- 
structed from the Green's and distribution functions at the middle of a closed 



4 



equilibrium SNS j unction . Furthermore, if d is much smaller than the coher- 
ence length ^0 = v^/i® /A, these Green's functions take the BCS form, with the 
phase-dependent proximity gap = A| cos((|)/2)| [18]. Within such model, 
the problem of current statistics in the Andreev interferometer reduces to the 
calculation of the CGF for an NS junction with the effective order parameter 
A^ in the superconducting reservoir. 

Proceeding to this calculation, we encounter a common difficulty, violation of 
the standard triangle form of G in the Keldysh space resulting from the gauge 
transformation in Eq. (5). In such situation, Eq. (4) cannot be decomposed 
into the Usadel equation for the Green's functions and the kinetic equation for 
the distribution functions, which makes well developed methods for solving 
Keldysh-Green's equations quite unusable. For this reason, the FCS problem 
in NS structures generally requires a numerical analysis of the whole 4x4 
matrix boundary problem which has been so far performed only in the limit of 
small characteristic energies {eV, T} <C A [9, 10]. 

In some particular cases, an analytical solution of this problem can be ob- 
tained by means of the generalized circuit theory [19, 20]. Within this approach, 
the CGF for a mesoscopic connector between two reservoirs is expressed in 
terms of the distribution p(T) of the transparencies of the conduction channels, 

Six) = ^JdEj^dTp{T)TT\nW{E,T,x), (6) 

W = 1 + (T/4)({Gz.,Gk(x)}-2), (7) 

where g is the connector conductivity. Equation (6) generally applies to the 
normally conducting structures with arbitrary p(T). It was also applied to the 
superconducting tunnel junctions [21] and point contacts [20, 22] with a singu- 
lar transparency distribution localized at the junction transparency. In general 
NS structures, due to dephasing between the electron and hole wavefunctions 
described by the commutator term in Eq. (4), statistics of conducting modes, 
in contrast to normal structures [23, 24], does not reduce to statistics of trans- 
parencies but concerns full scattering matrices. However, if the characteristic 
energies are much smaller than the Thouless energy, {eV, T} <C Ejh = h'D jlJ, 
the dephasing term in Eq. (4) can be neglected, which makes it possible to apply 
the transparency statistics for a normal wire [23] to the superconducting struc- 
ture. In long diffusive junctions, L » where the Thouless energy is small, 
Ejh <^ A, the quasiparticle spectrum is structureless at small energies, E <^ Ejh, 
which results in hnear voltage dependence of the CGF and, correspondingly, 
of all cumulants at eV <C Ejh [14]. In the opposite limit, eV S> Eju, the CGF 
for a long junction can be found within so-called "incoherent" approximation 
[25], by neglecting the contribution of the coherent proximity region. The cal- 
culations in [14, 25] lead to the conclusion that the FCS exhibits the reentrance 
effect: In both limits, eV <^ Ejh and eV » Ejh, it is described by the same ex- 
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pression for S{%). An interesting situation occurs in NS junctions with opaque 
interfaces dominating the net resistance of the junction [26]. In this case, the 
crossover between the coherent and incoherent transport regimes occurs at very 
small voltage of the order of the inverse dwell time of quasiparticles confined 
between the interface barriers. 

In this paper, we focus on the case of short NS junctions with the length 
smaller than and, correspondingly, with large Thouless energy. En ^ A^. 
In such situation, the energy region of negligibly small dephasing, E <C Ejh, 
overlaps with the region £ » A,],, in which the NS junction behaves as the normal 
system. This enables us to apply Eq. (6) and the transparency statistics for 
diffusive normal conductor at arbitrary voltages and temperatures, and present 
analytical solution of the problem. 

The calculation of the integrand in Eq. (6) is briefly discussed below. The 
Keldysh-Green's function Gr{x) in the normal reservoir is traceless in the 
Keldysh space and therefore it can be expanded over the Pauli matrices T as 

=^(gi+Ozgz), 8i8z = o, gi+gz = i, ^ = (i;;c,'i;y,Tz), (8) 

where the vectors gi^z{x) are expressed through local-equilibrium distribution 
functions in the voltage biased electrode. In the subgap energy region, £ < A^, 
the function Gl at the junction node is the unity matrix in the Keldysh space 
proportional to the Nambu matrix Green's function g in the superconductor, 

GL = | = avexp(a^0^), f = l, 0^ = arctanh(£/A^). (9) 

Then the calculation of the trace in the Nambu space in Eq. (7) is reduced to 
the summation over the eigenvalues a = it 1 of the matrix g, 

TrlnW = 5]^Tr^lnWa, Wa = a + zb, (10) 
a=l-T/2, b={T/2){cgi-ig,smhQ^). (11) 

Noticing that any 2x2 matrix can be presented in exponential form as 

Wa = exp{\nw + <pp), (12) 
w^ = a^ — W', cosh(p = a/w, p = xB/wsinhcp, Tr^ = 0, (13) 

where w is independent of a due to orthogonality of the vectors gi and g^, one 
easily obtains TrulnWo = Inw'^ and TrlnW = 21nw^. At £ > A,|,, the function 
Gl is traceless in the Keldysh space, 

GL = g{^gL), g = f^z^w{<^x%), e,^ = arctanh(A,|,/£), (14) 

where the vector gi is constructed from the equilibrium distribution function at 
zero potential. In this case, the 4 x 4 matrix W has the form W = a + db, where 
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a and b are scalars, 

a=l-(T/2)(l-gz.gzCOshe^), (15) 
P = (T/2)2 [(lifi )2 - (1^ X g,f sinh^ e<^] , (16) 

therefore it can also be transformed to the matrix exponent similar to Eqs. (12) 
and (13), with the traceless matrix p = dZ^/w sinh (p. Following this line, we 
obtain TrlnW = 2\nw^, and then, performing integration over T in Eq. (6), we 
arrive at the final expression for the CGF, 

where the quantities 9 and 9± are given by explicit relations, 

Z(0) cosh^ e = Z(2x) cosh^ 6^, (18) 
Z(0)coshe± = [Z(x)+cosX- l]coshe,^±tanh|[sinhp- (19) 



sinh {p - ix) - ismx]\ 1 - ^°^!^^^ , sinh^ 9,|„ 



9d, = arctanh 



cosh p—l 



(20) 



Z(x) = cosh(e) + cosh(p - ix), e = E/T, p = eV/T. (21) 

By using Eqs. (2) and (17)-(21), one can obtain analytical expressions for all 
cumulants. At zero temperature, the calculation essentially simphfies. Indeed, 
at r ^ and E > eV, the dominating terms in Eqs. (18)-(21) are proportional 
to exp(£), and therefore 9 and 9± are equal to Q^. This implies that the CGF 
is independent of the counting field at these energies, and all cumulants turn to 
zero. At £■ < eV, the terms with exp(;7 — inx) dominate, and we obtain 

cosh 9 = e~'^ cosh 9,^ , cosh 9± =e~'^ cosh Q^±{e~'^-l). (22) 

At subgap voltage, eV < A,^, when the charge transport at T = is only due 
to the Andreev reflection, the current /, the shot noise power Pj, and the third 
cumulant C3 read 

/ = lAq{z) , q{z)= T— arctanhx, Pi = 2e[l- /a/(z- ^ )] , (23) 

Jo X 

C3=N-^[{5z'-3)f{z-')+z], /A = ^, Fa = ^, (24) 
2z e e 

f{z) = (1/2) [z- (z^ - l)arctanhz-i], z = eV/A^. (25) 
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Figure 2. Shot noise power and third cumulant vs superconducting phase (a,b) at different 
voltages and J = 0, and vs voltage at different temperatures (c,d). Dashed lines denote voltage 
dependencies in the normal state at J = 0. In the panel (d), zero-bias slopes of the normalized 

C^iy) are indicated. 

At small voltage, eV <C A,|,, the shot noise doubles, P = (4/3)e/, and C3 = 
AN/\5 is four times larger compared to the normal case [8, 2, 27-1]. When 
the voltage increases and exceeds the gap edge, eV > A^, the normal electron 
processes at the energies E > A,^ begin to contribute to the charge transport, 
providing the normal-state voltage dependencies of the cumulants at eV » A^. 
At large voltage, the Andreev reflected particles produce voltage-independent 
excess components of the cumulants (in ballistic NS junctions, existence of the 
excess noise was first predicted in Ref. [28]), 

2 



1 = 



iN-hf{z)+r 
2 



2eh{z'-l)f{z)+Pr. lN = gV, 
(z + l){(z- l)[8z/3 - (8z^ - 3)/(z)] - 1/3} + cr, 



/a / 71 



-1 



Na 
2 



(26) 
(27) 

(28) 



At nonzero temperatures, T 7^ 0, we calculate the cumulant spectral densities 
/(£■), P{E) and C{E) defined as 

poo poo /•oo 

I = Ia dEI{E), Pi = 2elA / dEP{E), C3 = Fa / dEC{E). (29) 
Jo Jo Jo 
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Here /(£") = /isinh p/Z{0), and the functions P{E) and C{E) at £" < A(|, read 
P{E) = [2Qh + (1 - fi) sinhV] ,2 = 1 + coshecosh;?, (30) 



C(£) = 1^ [4/isinh2e + (2/2+3/3)sinh2;7+2e(3(l-/2)-2/i) 
whereas at £ > A,], they are given by equations, 
P{E) 



(31) 



6(1 + 2/1-2/2 '"';;^ , \ \+,M,^p-Z{0) 



Z2(0) L^^V' ' ^"'coshs+l ' ' ^ ^^"^ ' ^^^^ 

^(^^ = z3(0)(i+coshe) {"^^^^+'"'^^)(g + '"^'^^+ (^^^ 
3 [z(0) ( 1 - 2/3 ) + e (4 ( 1 - /2 + /3 cosh s) + 3 cosh 8 - 2/3 ) + 

sinh^S^2/3 — coshs+ (3 — 5coshp)/2^ +/2(5coshs— l)sinh^/7| , 

/i =e^cothe^, /2 = (/i- l)/sinh2e<^, /s = (/2- l/3)/sinh2e^. (34) 

In Eqs. (30)-(34), the functions /■ describe energy variation of quasiparticle 
spectrum which is most essential in the vicinity of the gap edge A(|,. 

As shown in Fig. 2,(a,b), the cumulants oscillate with the phase and exhibit 
deep minima at ^ mod 27i = 71, when the gap closes and the cumulants approach 
their normal values. When the proximity gap A,], approaches eV , Pi{^) exhibits 
a peak, while Ct,{^) shows a step-Uke structure. The voltage dependence of 
the cumulants for different temperatures is plotted in Fig. 2,(c,d) in specifically 
normahzed variables, which provides universality of the curves for any (j). As 
the temperature increases, the current noise approaches finite value at = 
due to thermal fluctuations, and exhibits quadratic dependence on the applied 
voltage at eV <^ T. Within the intermediate voltage region, T <eV < A^, Pi{V) 
becomes linear with doubled slope produced by the Andreev reflected particles, 
and at eV > A^, the slope turns to its normal-metal value. A considerable 
excess noise at large voltages is contributed by both the thermal fluctuations 
and Andreev reflection. A more interesting behavior is discovered for the third 
cumulant. As the temperature departs from zero, the zero-bias slope of the 
normahzed Cj{V) increases by factor of 5 compared to the zero temperature 
(which is similar to the normal structure [4]), and approaches the value 4/3. 
Then, atT <eV < A,^, the slope of the normalized C3(V) returns to the subgap 
value 4/15 for r = 0. AteV ^ A,^, the curve C3{V) shows 5V!-like feature, and 
finally, at eV > A,^, it approaches a straight hne with the (normal-state) slope 
1/15. This implies that C3 acquires anomalously large thermal component at 
voltage eV ^ Ai^, which, however, rapidly decreases at eV > A,^ and/or T > A,|, 
towards the normal metal level. 
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Figure 3. Effective transferred charge (a), differential noise (b), and differential resistance 
(inset) vs voltage at (|) = (solid lines) and (|) = 0.731 (dashed Unes), T=0. 

The voltage dependence of the effective transferred charge defined as qeff = 
{3/2)dPi /dl [10] and the normalized differential noise (3/?/2) dPj/dV is shown 
in Fig. 3 for r = 0. At zero voltage, both quantities approach 2e, whereas at 
eV » they turn to e. However, the behavior of these quantities becomes 
very different when the voltage approaches the proximity gap edge A(|,: the 
normalized differential noise increases to 3e, whereas the effective charge q^ff ~ 
{dP]/dV){dV /dl) rapidly turns to zero, along with the differential resistance 
dV/dl (see inset in Fig. 3). We note that the effect of zero differential resis- 
tance giving qeff = results from the resonant transparency enhancement at 
the proximity gap edge, due to the singularity in the density of states. Such an 
"exotic" behavior of these quantities at the gap edge leads us to the conclusion 
that none of them can be unambiguously associated with the physical elementary 
transferred charge, in contrast to what is commonly suggested. Similar effects 
have been predicted in Ref. [26] for an NS structure with opaque interfaces; 
however, in such case, a considerable enhancement of dPj/dV and suppression 
of dPj/dl occur at small applied voltage of the order of the inverse dwell time 
of quasiparticles. 

It is instructive to compare our analytical results for short-arm interferometers 
with that obtained numerically for long NS junctions with a small minigap 
Ef, ~ Exh ^ A. The results are qualitatively similar: in long junctions, q^ff' is 
equal to 2e at eV <C Eg and has a minimum at eV « Eg{i^), which moves towards 
small voltage when at (|) — >^ 7t [10]; the differential noise is also non-monotonous 
and approaches maximum at eV ^ SEjh [9]. After this comparison we see that 
the proximity gap A^ in short junctions plays the role of the minigap Eg in long 
junctions and determines the feature in the effective charge, though this feature 
at eV ~ Eg in long junctions is much less pronounced. However, as noted above, 
a qualitative difference of long junctions is the existence of an intermediate 
incoherent voltage region Eg <^ eV <^ A, where both the effective charge and 
the normalized differential noise have the value 2e, and their crossover to e 
occurs only ateV >A [29]. 
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Appendix 

For reference purposes, in this Appendix we outline the procedure and summarize the results 
of calculation of the CGF for a diffusive connector between normal reservoirs, by using the 
extended Keldysh-Green's technique. For generality, we consider a diffusive wire interrupted 
by tunnel barriers, which enables us to present several original results and to examine various 
limiting situations. 

In normal systems, the matrices G and / are traceless in the Keldysh space and therefore they 
can be expressed through 3-vectors with the components diagonal in the Nambu space, G = gz, 
I = ri, where x is the vector of the matrices x, and = I. Since the left-hand side of Eq. (4) 
turns to zero in normal systems, the formal solution of Eq. (4) for the matrix current density 
in each segment of the wire can be easily obtained, 

4f = givlnGiG2=^ivln[gig2 + !^(giXg2)]=^4f, lN = igNP<^N, (A.l) 

where gj^ is the conductivity of the wire segment, G12 are the Green's functions at the left and 
right segment edges, respectively, (^ff = arccos^ig2 is the angle between the (complex) unit 
vectors gi and g2, and p = {gi x gj()/ s'm(^ff is the unit vector perpendicular to ^1 and g2- 

The matrix current Is through the tunnel barrier can be expressed through the Green's func- 
tions G- and G+ at the left and right sides of the barrier by using the boundary condition [30], 

/b_ =/b+ = y[G_,G+] =x/ij, lB = igBPsin(^B, <|)b = arccosg_g+, (A.2) 

where p = {g^ x ^+)/ sin(|)g and gB is the barrier conductivity. 

The conservation of the matrix current along the connector, / = const, following from Eq. (4) 
and the boundary condition in Eq. (A.2), results in conservation of the vector current, 7 = /jv = 
Ib = const. This implies that for all elements of the connector, the unit vectors p coincide, 
therefore the Green's vectors g lye within the same plane and the vector p can be constructed 
from known Green's vectors gi and gR in the reservoirs, p = {gL x gij)/sin(|), where (j) is the 
angle between g^ and gR. From the current conservation, we also conclude that all barriers and 
wires are characterized by a single variable r\, 

gB sin(|)B = gffitiN = gil = const, (A.3) 

where the normalization constant g is chosen to be equal to the conductance of the whole 
connector. Thus, the vector current is given by equation, 

I=^{gLXgR)- (A.4) 
sm(|) 

The planar rotation of the Green's vector along the connector results in the additivity of the 
angles between all consecutive vectors g, therefore the sum of these angles is equal to (|), 

52 'I'A' + E '^B = '^ = atccos{gigR), (A.5) 
wires barriers 

Mr = l+P-+{e^-l)+ P+- (e-"' - 1) , (A.6) 

-Poa' =«c(l-«a'); n-=nF{E), n+=nF{E + eV), (A.7) 
which leads to the following equation for the parameter T|((|)), 
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where Rn is the net resistance of all wires, Rk is the resistance of the k-ih barrier, and R = g . 
By using the definitions in Eq. (3), we obtain the counting electric current /(%) and the CGF, 

I{X) = ^ TdE TtoJ, = r dE Tr ^ {gL xgR),, (A.9) 
2e Jq 2e Jq sm(j) 

Six) = ^jdETr [r^li2/2 + 1, ( ^ " 4^^^) ' 

We note that the statistics is insensitive to the position of the barriers, and depends only on the 
barrier resistances and the net resistance of the diffusive part of the connector. In the absence of 
barriers, rj^ — > 0, the CGF has the form. 

Six) = / dEi^'^ = ^jdE arccos^ [l {^S - l) +P+- [e''^ - l)] . (A.ll) 

At zero temperature, the integration over energy in Eq. (A. 10) can be explicitly performed, 

Six) = I [rW/2 + i:, (l-i/l-r2Ti2)/r,], (A.12) 
where N = gVto / e. From Eq. (A. 12) we find the Fano factor F in the shot noise power P/ = eFI, 
F = (2/3)(1 + 2S3), S„=£^,f, (A.13) 

which varies between the Poissonian value F = 2 for the tunnel connector and 1 /3-suppressed 
value, F = 2/3, in the absence of barriers. The third cumulant C3 varies between A' for Poissonian 
statistics in the single barrier case and Af/15 for a diffusive conductor, 

C3(y,0) = (AT/IS) [1 + IOS3 (1 +4S3) - 36fi5] • (A.14) 

It is interesting to note that Eq. (A.8) can be easily transformed into equation for the trans- 
parency distribution p (T) , by making use of the relation of the generalized circuit theory between 
the counting current /(x) and the matrix current / following from Eqs. (6) and (3), 

'^'^^ = Te[ ^^=\f^dlpil)i:\GL,GRix)]W-'. (A.15) 

Rewriting this equation in the vector representation, comparing it with Eq. (A.9), and introducing 
the variable z = (1/2)(1 — gigR), we obtain the integral equation for p(T), 

/'I^ = _J=, (A.16) 
io 1-zT 2y^(T^' 

where t] obeys Eq. (A.8) with the function (|) = 2arcsin^ in the right-hand side (rhs). The 
solution of Eq. (A.16) has the form p(T) = Rer|/27rTVl -T, where Ti(T) is the solution of 
Eq. (A.8) with the function 7t + 2(arccosh(l/\/T) in the rhs [24]. 

In some limiting cases, one can obtain an analytical solution of Eq. (A.8). In particular, if the 
number M of the barriers is large, M S> 1, then the resistance of each barrier is small compared 
to the net resistance, Rk -C R. In this case, the approximate solution of Eq. (A.8) is T) = (|), and 
the CGF coincides with that for diffusive wire, 5(x) =N arccos^ e''^/^ . In the tunnel limit, when 
the resistance of each barrier much exceeds the net resistance of diffusive segments, R^ ^ Rn, 
the first term in Eq. (A.8) can be neglected. Then an analytical expression for the parameter T] 
and the CGF at arbitrary M can be obtained in the case of equivalent barriers, rk = \ /M, 
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when the Fano factor is given by F = (2/3) (l +2/M^). In the Umit of large number of the 
barriers, 1, we return to the "diffusive" statistics, while for single-barrier structure, M = 1, 
we obtain pure Poissonian statistics, 5(x) = N{e''^ — !)• 

At arbitrary temperature, the cumulants can be found analytically by asymptotic expansion 
in Eqs. (A.8), (A. 10) over small T| and %. In particular, the noise power, 

i'/(V,r) = ^ [(l+2S3)|coth|+2(l-S3)] , (A.18) 

exhibits crossover between the shot noise at r ^ eV and the Johnson thermal noise Pj = 4r//? 
at large temperature, T » eV. The voltage dependence of the third cumulant, 

C3 (F, r) = C3 ( V, 0) + |iv( 1 - lofi^ + 9fi5 ) '"t!frL ' 

5 psinh^(p/2) 

is linear in both limits and approaches (A'/3)(l +263) at high temperatures. In the limit of 
tunnel connector, S„ = 1, the second term in Eq. (A. 19) vanishes, and C3 becomes temperature 
independent. In the absence of barriers, S„ = 0, Eq. (A. 19) reproduces the result of a modified 
kinetic theory of fluctuations for a diffusive wire[4]. 

In order to access PCS in multi-terminal structures, which consist of a set of connectors 
attached between several normal electrodes and a diffusive island (node) with negligibly small 
resist£ince, separate counting fields Xa and parameters T|a are to be introduced in each arm[16], 

la = ikaiSa, X gc), ka = gaTla/sin(|)a. (A.20) 

The quantities y\a obey the equations similar to Eq. (A.8), with the angles (|)a = arccos(^cift ) 
in the rhs, where the Green's vector gc at the node can be found from the current conservation 
law, £„ 4 = 0, 

g, = G|^^&■, G = Y,^^aga- (A.21) 

According to Eq. (A.21), the vector gc depends on all counting fields Xa. which reflects cross- 
correlations between the currents in different connectors. Por the system of tunnel connectors, 
where the quantities i,a are equal to the conductances ga and therefore become independent of 
X, the CGP at zero temperature can be explicitly evaluated [21], 

S{X} = ^G^l + Al^^gyg^ie^Xa - 1), g„ = g„/G, G = £pgp, (A.22) 

where the index V denotes the voltage biased electrode. 

Por arbitrary connectors, the cumulants can be found from asymptotic solutions of the equa- 
tions for Tia and gc at small Xa- For instance, the partial current through a-th connector is 
la = Vgagy, and the Fano factors defined as F^p = (2ei'//a)(3/a{x}/9X|3)x=o read 

^ap = (2 - hv) Sap - ^Ip [1 + Iv (S3a + Ssp) - B3v{i -gv? -gv L Sy^^y] ' (^-23) 

The diagonal elements Faa of the matrix F^p have the meaning of the Fano factors for the 
shot noise in a-th connector and may vary between 2/3 and 2, approaching Poissonian value 
^aa = 2 for a large munber of the terminals, when the normalized conductances g„ become 
small. The cross-correlators F^p (a P) between the currents in different terminals are negative 
[29]. In a particular case of diffusive connectors (Ba = 0), Eq. (A.23) reproduces the result of 
Ref. [31] for a so-called star-shaped geometry. 
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